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We present a theory for the temperature-dependent nuclear magnetic shielding ten- 
sor of molecules with arbitrary electronic structure. The theory is a generalization 
of Ramsey's theory for closed-shell molecules. The shielding tensor is defined as a 
second derivative of the Helmholtz free energy of the electron system in equilibrium 
with the applied magnetic field and the nuclear magnetic moments. This derivative is 

P-i 

analytically evaluated and expressed as a sum over states formula. Special consider- 
ation is given to a system with an isolated degenerate ground state for which the size 
of the degeneracy and the composition of the wave functions are arbitrary. In this 
case the paramagnetic part of the shielding tensor is expressed in terms of the g and 
A tensors of the EPR spin Hamiltonian of the degenerate state. As an illustration 
^ of the proposed theory, we provide an explicit formula for the paramagnetic shift of 

the central lanthanide ion in endofullerenes Ln@Cgo, with Ln=Ce 3+ , Nd 3+ , Sm 3+ , 
Dy 3+ , Er 3+ and Yb 3+ , where the ground state can be a strongly spin-orbit coupled 
icosahedral sextet for which the paramagnetic shift cannot be described by previous 
theories. 
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I. INTRODUCTION 



Nuclear magnetic shielding measured in NMR spectroscopy is caused by the electrons' 
response to an externally applied magnetic field. In closed-shell molecules this response 
consists of an induced magnetic field caused by electronic orbital currents, the electron spin 
playing no role other than forcing appropriate orbital permutation symmetries in the many- 
body wave function. This, at least, is true when spin-orbit coupling is ignored. Ramsey, 
in Ref. Q, presented the quantum mechanical theory of the shielding tensor for molecules in 
this regime. 

In open-shell molecules however, there is a second contribution to the shielding, arising 
from the permanent magnetic moment associated with the spin of the unpaired electrons. 
This contribution is known as the 'paramagnetic shift'.- - — More generally, a paramagnetic 
shift arises from a degenerate electronic state. When the degeneracy is weakly split in 
a magnetic field, the magnetic polarization of each thermally populated state can be de- 
scribed in terms of permanent non-compensating spin (and orbital) currents,- inducing a 
net response field at the nucleus, each state contributing in proportion to its thermal popu- 
lation. A general equation for the isotropic paramagnetic shift was first derived by Kurland 
and McGarvey.- Moon and Patchkovskii derived a formula for the paramagnetic shielding 
tensor in an arbitrary Kramers doublet and expressed this formula in terms of the g and 
A tensors of the EPR spin Hamiltonian. 5 Later Pennanen and Vaara extended the Moon 
and Patchkovskii theory of paramagnetic shifts to cases which deviate from a pure spin 
degeneracy in the lowest order of perturbation theory in the spin-orbit coupling. 7 

Although the work of Pennanen and Vaara represents the first comprehensive effort to- 
wards a completely general theory of NMR chemical shielding for electronically degenerate 
states in the weak spin-orbit coupling regime, their approach- exposes a few issues: (i) It 
is developed using perturbation theory on the thermodynamic internal energy U, which 
appears to require different perturbation expansions of a same energy level, according to 
whether the level is either a weighted addend, or an exponent of the thermal weight, in 
the Boltzmann sum-over-states, (ii) It considers Zeeman and hyperfine spin Hamiltonians 
linear in the spin operators, which is strictly correct only if the degeneracy is no larger than 
threefold, (iii) It does not provide a formulation of the problem for an arbitrary electronic 
degeneracy, such as can be found for example in strongly spin-orbit coupled systems, (iv) 
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Finally, the inclusion of low-lying excited states resulting from weak splitting of a degener- 
acy, as described by a zero field splitting Hamiltonian for pure spin states, is not described 
in the most general way. 

We have recently published a communication in which especially points (iii) and (iv) were 
thoroughly discussed and generalized.— In the present paper we wish to present a theory 
which offers a rigorous solution to point (i). We do this by taking no more than Ramsey's 
original assumptions^ and applying them to a system where the electrons are in thermal 
equilibrium. Thus we are led to identifying the shielding tensor with a second derivative of 
the Helmholtz free energy, as in Eq. ([3]). Then, without needing further assumptions, we 
use perturbation theory on the free energy to derive the expression for the shielding tensor, 
Eq. (fTOl) . In the second part of the paper, we show that this expression lends itself easily 
to a reformulation of the paramagnetic shielding in a degenerate state in terms of the spin 
Hamiltonian g and A parameters (Eq. f)16p ). which provides a general solution to point (ii). 

II. NON-DEGENERATE GROUND STATE: SHIELDINGS AS ENERGY 
DERIVATIVES 

In the following we shall be concerned exclusively with nuclear shielding in the 'solid 
state limit',- by which is meant that the molecular nuclei are fixed in space with respect to 
the static external magnetic field. 

In this section we want to recall the essential points of Ramsey's theory for a molecule in 
a non-degenerate ground state.- He observed that the nuclear magnetic moment fi can be 
treated as a classical vector because it is much slower in its dynamics than the electrons. The 
idea is that as the nuclear moments go about their rotating motion the electron cloud follows 
adiabatically; at every moment the electrons are in the ground state that corresponds to 
the instantaneous orientation of the nuclear moments. Hence the components of /j, become 
external parameters in the electronic Born-Oppenheimer Hamiltonian, as are the nuclear 
positions and the applied field B. The ground state energy will naturally be a function of 
/x and B, and can be expanded in a Taylor series as 



Here E is the electronic energy in the absence of nuclear moments and external field. Note 




(1) 
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that the energy contains no terms of odd degree, which is a consequence of the assumption 

that Eq is a non-degenerate eigenvalue of a time-even Hamiltonian. The second term in ([1]) 

has the form of a Zeeman interaction between the nuclear moment fi and the induced field 

— B ■ (7. The 3x3 matrix cr is the magnetic shielding tensor of the nucleus. We consider 

a to be field independent, but note that field-dependent corrections to the shielding can be 

found among the higher order terms of (pQ). It follows from Eq. ([1]) that can be obtained 

from perturbation theory (up to second order) on Eq. Thus, in Ramsey's theory, 

3 2 E 
%3 dBid/jj 

For reasons of simplicity we have assumed that there is only one magnetic nucleus in the 
molecule. When there is more than one, Eq. (CTJ is readily adjusted to expand the energy 
E(nx, . . . , B) in all the nuclear moments. Each nucleus will have its own shielding tensor. 

III. ARBITRARY ELECTRONIC SPECTRUM: SHIELDINGS AS FREE 
ENERGY DERIVATIVES 

A. Separation of slow and fast dynamical variables: Free energy as an 
effective Hamiltonian for nuclear spins 

We have seen that in the case of a non-degenerate isolated ground state the shielding 
tensor is obtained from the electronic energy. In this section we want to extend the theory 
to include molecules which have several thermally populated electronic energy levels, so that 
temperature comes into play. We want to find an expression for the shielding tensor that 
is universally valid, regardless of the particular electronic spectrum of the molecule. We 
start again from the physical assumptions of Ramsey: the nuclear magnetic moments are 
classical vectors, they interact with the electrons, and the electron cloud is at all instants 
in a state of equilibrium governed by the applied magnetic field and the instantaneous 
orientation of the nuclear moments. This equilibrium state is not simply the ground state, 
because we have to allow for the molecule to possess multiple electronic states between 
which rapid transitions occur. As the timescale of these transitions is much shorter than 
the timescale of the nuclear spin dynamics,- 1 ^ the electron system can at all times preserve 
thermal equilibrium at constant temperature T = 1/(3. That is, the state of the system is 
described by a Boltzmann distribution over the energy levels E n (/j,, B), eigenvalues of the 
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electronic Hamiltonian i?(/x,B). As the energy levels depend on /j,, so do the partition 
function and all other thermodynamic functions. Note that /j, is not treated as a dynamic 
variable of the thermodynamic system (which consists of the electrons only) but as an 
external parameter that modifies the energy levels of the system. 

Our next task is to establish which equation should replace Eq. (Q. Recall that Eq. ([1]) 
gives the electronic energy as a function of nuclear moments and external field. For the 
present purpose however, a better interpretation is to regard this energy as the effective 
Hamiltonian 11 for the system of nuclear moments in the external field B. There is no longer 
any explicit reference to the electrons but the influence of the electrons on the motion of the 
nuclear moments enters the effective Hamiltonian in such quantities as the shielding tensor, 
the nuclear spin-spin coupling, etc. For the class of molecules to which Eq. (JTJ) applies (non- 
degenerate, isolated ground state) we can thus write: H e g(fj,; B) = E(fjt, B). Note that in the 
effective Hamiltonian, /x is a dynamical variable, whereas B remains an external parameter 
(hence the semicolon). In the more general case of electrons in thermal equilibrium the 
effective Hamiltonian is temperature-dependent and is givenii by the electronic Helmholtz 
free energy F = U — TS, where F, U, and S are parametric functions of n and B, and of 
temperature: 

if e fj(/x; B) = F(fjt, B) = F + BiCTijfij + higher order terms. (2) 

ij 

Here we have expanded F in a Taylor series in analogy with the expansion of the ground state 
energy in Eq. ([T]). One way to see that F acts as the effective Hamiltonian for the nuclear 
moments is by considering the work done on the system by changing an external parameter 
(/ij, say) at constant temperature. This amount of work is given by the concomitant negative 
change in H c g. The only thermodynamic function that has this property is the Helmholtz 
free energy. Hence, up to an irrelevant constant, the effective Hamiltonian is given by F. 

B. Analytical derivatives of the Helmholtz free energy 

The shielding tensor component is the coefficient of Bifij in the effective nuclear 
Hamiltonian. In Ramsey's theory is temperature independent and can be obtained 
from Rayleigh-Schrodinger perturbation theory on the ground state energy, Eq. ([1]). In the 
general theory Cy is temperature dependent and can be obtained from perturbation theory 
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on the free energy. From Eq. (|2j) we have 

d 2 F 



(3) 
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We write the electronic Hamiltonian B) = Hq + V(fi, B). Here Ho is the electronic 
Born-Oppenheimer Hamiltonian in the absence of nuclear magnetic moments and external 
magnetic field, and V(fi, B) collects all interactions of the electrons with \i and B. The 
expressions of these terms are well known and can be found for example in Abragam.— 
Their detailed form is not important for our discussion. 

In the calculation of o^-, only those terms of V will contribute that are linear in or 
linear in B or bilinear in /x and B. Thus we write V as the sum of four parts: 

V z = -m ■ B = - 2J mBi, 

i 

Vhf = T ■ A* = Tilh, 

(4) 

ij 

V' = terms of higher order in B and fi. 

Here V z is the electronic Zeeman Hamiltonian, with m = — (L + g e S); Vm is the hyperfine 
coupling, and J- may be further divided into orbit, spin dipole, and Fermi contact terms;— 
Vp is the diamagnetic nuclear-magnetic field coupling; V contains all those terms that are 
not needed to calculate the shielding tensor, and will be discarded from now on. In general, 
Vhf and Vr> must be summed over all magnetic nuclei, each with its own J- and T>, as these 
depend on the position of the nucleus which they represent. It is important to keep in mind 
that m, and T> are electron operators, whereas and B are classical external parameters. 
Furthermore, m and J- are time-odd, whereas T> is time-even. 

For the derivation 13 of a perturbation expansion of F it will be convenient to introduce 
a parameter A in the Hamiltonian to denote the combined order in fj, and B: 

H — H + XVx + AV 2 , with Vi = V z + V m and V 2 = V v (5) 

The free energy is given by 

F = — — InTr p, where p = e~^ H 
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Now p can be expressed as a power series in A: 

p = p + Api + \ 2 p 2 + 



(6) 



Here po = e an d Pi and p2 are first and second order corrections. A power series for F 
can now be obtained as follows: 

F — — — ln(Trp + A Tr pi + A 2 Trp 2 + . . .) 







1 , m 1 , / > Tr p-i 

--lnTrp --ln 1 + A-^i + A^ 
/3 /3 V Tr Po 



A Tr pi A 2 
PT^p~ + J 



1 /Trp] 



2 VTrpo 



Trp 2 
Trpo 

Trp 2 
Trpo 



+ 



+ . . . . 



We know however that the free energy is an even function of the time-odd fields B and /x 
and therefore odd powers of A must not occur in the expansion of F. Hence Tr p\ = and 

A 2 Trp 2 



Trp 



+ 



(7) 



To proceed, an expression for p 2 is needed, p satisfies the Schrodinger equation || = —Hp. 
By substituting Eqs. ([5]) and ([6]) and collecting terms of equal power in A, a recurrence 
relation for pi is found: 

dpi 

d(3 



-H p 1 - Vipo, 



^ = -H p 2 - Vipi - V 2 p , 

etc. These are linear differential equations of the form ^| = ay+b(x) } whose general solutions 
are given by y = e ax e~ at b{t)dt. Thus we find 

pi = - / e {w -P )H »V ie - wHo dw 
Jo 

p 2 = - e ( - w -P )Ho V 2 e- wHo dw+ / e^-^V^'-^V^-^dw' dw. 
Jo Jo Jo 

Taking the trace of p 2 , we use Tr(AB) = Tr(BA) to rearrange and simplify the integrands. 13 



Write p 2 = pp + pp , then 



Trp 2 



(i) 



where 



Tr[e-P Ho V 2 )dw = -P(V 2 } Trp , 
Tr(r 2 po) 



(V 2 )o 



Trpo 



is the thermal average of V<i in the canonical ensemble corresponding to Hq. For the second 
part of p2 we have 

Trpj, 2) = / / Ti[e-^ Ho e^-^^V^ 7 "' -^V^dw' dw. 
Jo Jo 

The integrand depends on w — w' only, suggesting a change of variables u = w, v = w — w'. 
This gives 

Trp®=[ [ Tile-^e^Vxe-^V^dvdu. 
Jo Jo 

A second change of variables u — > ft — u, v — > ft — v shows that this integral is also equal to 
Jo I °^ ^ ne same integrand. Hence 

Trp { 2 ] = - / Tr[e- f>Ho e vHo V 1 e- vHo Vi]dvdu 
2 Jo Jo 

= 2 (/ e vHo V 1 e- vB °Vidv^ Tr Po . 
At this point A is no longer needed and we set A = 1, so that Eq. ([7]) becomes 
F = F + (V 2 ) - \ Qy^V.e-^V.dw^ +.... 

We substitute V\ and V 2 according to Eqs. ([5]) and @ and then take the second derivative 
as in Eq. ([3]) to obtain the shielding 
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(A;)o + ( / e^rme-^^dw ) . (8) 
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In deriving this expression we have used that, for any two operators A and B, 

(/ ] eWH ° Ae ~ WH ° Bdw ') = (J Q ' e wHo Be- wHo Adw^ , 

which can be shown by a change of integration variable w — >■ ft — w. 

To evaluate the integral in Eq. (|8]), the ensemble averaging is carried out before the 
integration. Let |nz/) be the eigenstates of H with eigenvalue E n (0,0), hereafter simply 
written E n . The index v labels an arbitrary orthonormal basis of states with the same 
energy E n . The integral then becomes 

rP 

/ (e wHo m i e- v ' ao J 7 j ) dw = 
Jo 

1 f 13 
— V, (nv\m i \mp)(mp\J r j \ni')e~ pEn / e w{En ~ Em) dw, 
Qo Jo 

nv,mfj, u 
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where Qo = Tr p Q = J2nu e ^ En denotes the partition function of H . The remaining integral 
can now be evaluated and this gives 



rP i r 

/ (e wHo m i e- wHo T j ) dw= — Ve* 0y2(nu\m i \nu')(nu'\J : j \nu) 

(n i/|mi|m n)(m y^T^n v) + c.c. 

1? _ 7? 



(9) 



Finally, combining Eqs. (jHJ) and ([9]), we find a sum over states formula for the nuclear 
shielding tensor, 



<Jij = — e /3 v') (n v\ J-j\n v) + ^~\(n v\T>ij\n v) 

(n v\mi\m /j,)(m n\J : j\n v) + c.c. 

/ J - / m J -'n 

m=f=n u,fi 



(10) 



Incidentally, the complex conjugate (c.c.) is not strictly needed and may be replaced by a 
factor of 2 in front of the last sum, because this sum is always real. The present expression 
bears an evident similarity to the Van Vleck equation for the magnetic susceptibility. 14 One 
only has to substitute rrij for Tj and the diamagnetic susceptibility term for T>ij to obtain the 
Van Vleck equation. The origin of this similarity is readily understood when one notices that 
the susceptibility tensor, just like the shielding tensor, occurs in a term of the free energy 
expansion Eq. (j2J), viz. in the term that is quadratic in the external field: BiXijBj. 
The calculation of Xij then proceeds along exactly identical lines and yields the Van Vleck 
equation. 

Eq. ( ITUl) can be cast in a more transparent form using projection operators. Let P n = 
s £2iJ\nv)(nv\ be the projector on level n (which is possibly degenerate), and Q n = 1 — P n 
its complement, then we can write 



fjP n TTl{P n J~^ jPn Pri^ij Pn 



Qn t~ r> r> ~r Qn 



~\~ PnT^i jt t-t 3~ jPn ~\~ Pn*Pj -rj j-, TTliPn 



This formula shows that the shielding tensor can be seen as a Boltzmann average over the 
electronic energy levels of the unperturbed molecule, of a quantity (the trace of the operator 
in brackets) that is associated with each level. This quantity has a temperature-dependent 



part and a temperature-independent part. The temperature-dependent part, 



(3Ti(P nmi P n ^P n ), 

which we shall call the Curie term, can only be non-zero if E n is degenerate. This is because 
if E n is not degenerate, the eigenstate \n) is real (w.r.t. time-reversal) and expectation 
values of time-odd operators on a real state vanish (i.e, P n rriiP ri = P n TjP n = if E n is 
not degenerate). Physically, the Curie term arises from the polarization in the applied field 
of the permanent magnetic moment of the degenerate level. The temperature-independent 
part of the quantity in brackets consists of two terms, the 'diamagnetic' term (~ T>) and 
the 'paramagnetic' term. These terms are known from the Ramsey theory for shielding in a 
non-degenerate ground state. When the state is degenerate, the Curie term appears on top 
of the Ramsey terms, shifting the chemical shift from the value it would have if the state 
were not degenerate. This shift, coming from the Curie term, is known as the paramagnetic 
shift. 

As Eq. ([TO]) was derived for a system with an arbitrary electronic spectrum, it must 
obviously also be valid for a 'closed shell' system, i.e., for a system with a non-degenerate, 
isolated ground state. In that case, Eq. ( TTTJT) reduces to 

^ww^^y^ (id 

nu n 

which is indeed the Ramsey expression, which can also be obtained starting from Eq. ([1]) 
and applying perturbation theory (first order in X> and second order in 7- and m) on the 
ground state energy. The equivalence of the two approaches, one starting from the ground 
state energy E, the other from the free energy F, follows from the fact that the entropy is 
zero (as only the ground state is occupied), and therefore F = E. We would like to note that 
although Eq. (fTTj) has the same form as the original Ramsey expression,- it is more general 
than the latter in that the ground state |0) in Ramsey's paper was assumed to be a pure 
spin singlet state, which allowed him to ignore the spin-dependent parts of J- and m. The 
present treatment imposes no restrictions on |0) other than those stated before, viz. that 
|0) is the non-degenerate ground state of a time-even Hamiltonian Hq (which, incidentally, 
excludes odd electron systems). H may include spin-orbit coupling and in cases where 
this is an important effect the spin-dependent terms, the spin Zeeman term in m and the 
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spin-dipolar and Fermi contact hyperfine terms in must be included in the calculation of 
<7ij from Eq. ( ITT]) . 

It should be noted that the procedure outlined above for the shielding tensor can be 
used to obtain expressions for any desired term in the nuclear effective Hamiltonian. It 
is, for example, straightforward to obtain nuclear spin-spin coupling tensors by using the 
appropriate V terms in the Hamiltonian and carrying out the subsequent derivation as 
before. Thus, in analogy with Eq. (jSJ), the coupling tensor between two nuclei K and L is 
found to be 



d 2 F 

J« L - 



3 

where T>^ L is the diamagnetic nucleus-nucleus coupling term.— 

IV. SHIELDING TENSOR IN A DEGENERATE ISOLATED GROUND 
STATE 

We now consider a case of particular interest: a system with a degenerate ground state 
of multiplicity u. We assume that the system is in thermal equilibrium at a temperature 



such that no excitec 
published in Ref. [10 



states are thermally populated. The main results of this section were 
, Here and in the appendices we provide a thorough discussion of the 
details of their derivation. 

Applying our general equation (ITU]) to the case of a degenerate and isolated ground state 
we find 

1 



CTij 



^^(OHmilOz/'XOz/'lJ^IOz/) 

u,u' 

V m^O v,n m " 

The last two terms are familiar from the Ramsey expression Eq. ( TlTi) (and reduce to the 
latter when oj — 1), but here they are averaged over the states of the manifold. The first 
term is the Curie term which, as we have noted before, is unique to degenerate states and 
causes a paramagnetic shift, inversely proportional to temperature. As much of the following 
discussion will focus on this term we repeat it explicitly: 

a% = ^^(0i/|7ni|0i/)<0i/|^|0i/). (12) 
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The sum in Eq. (|T2|) is evidently the trace of the product of two matrices, respectively the 
representations of the electronic magnetic moment and the hyperfme field on the nucleus in 
a basis of the degenerate manifold. Being a trace of a matrix product, it is invariant under 
unitary transformations of the basis, which is of course required for an observable property. 
This means that we are free to choose whichever orthonormal basis we want to calculate er?- 
using Eq. ([121). 

Both matrices in Eq. ffT2|) are known from the theory of electron paramagnetic resonance 
(EPR) spectroscopy. They represent respectively the Zeeman effect and the hyperfme cou- 
pling in the degenerate manifold. And thus we arrive at the well-known result that the 
paramagnetic shielding tensor can be obtained completely from EPR parameters, and, vice 
versa, NMR experiments on paramagnetic molecules provide information on the EPR pa- 
rameters. The precise expression of this correspondence will be derived in section HVBt but 
let us start with an example, the simplest case of degeneracy, a Kramers doublet. 



A. Kramers doublet 



In a Kramers doublet the Zeeman and hyperfine interactions can be conveniently ex- 
pressed by means of a spin Hamiltonian 

# 5 = /iBS-g-B + S-A-I (13) 

operating within a fictitious spin-1/2 doublet.— What this means is that the matrix elements 
of V z + Vhf (Eq. (jl])) within the true wave functions of the Kramers doublet are the same as 
the matrix elements of Hs within the fictitious spin doublet.— 1 ^ For example, if we denote 
the basis states of the Kramers doublet by |0 a) and |0 6) and let |0 a) correspond to the 
fictitious spin state ||) and |0 6) to | — |), then we have 

(0a\m x \0b) = -^ B (f Eiffel ~ \) = ~^{9xx ~ «flw)/2, 
(OapgOfc) = —(l^SiA^ -\) = — (A xx -iA yx )/2, 

and so on. Here gj is the nuclear ^-factor and we have used /x = gifi^l, where I is the 
nuclear spin vector in units of h. Using the spin Hamiltonian of the Kramers doublet, Eq. 
(fT2l) becomes 

a l = -7r^—J2 9kiAl i Tl ( SkSl ^ 
12 



The trace of the product of two cartesian spin components is given by 

And with S = 1/2, we find that the paramagnetic shielding tensor for a Kramers doublet 
can be expressed as follows. 

- P = -^g T A. (14) 

This equation was first obtained by Moon and Patchkovskii using a different method.- Note 
how the two interactions that combine to give rise to a paramagnetic shift are separated 
in Eq. (|T4"|) : g represents the Zeeman splitting in the external field and is a property of 
the molecule as a whole, independent of the nuclear moments; A represents the hyperfine 
interaction between the electrons and the nucleus whose shielding tensor we wish to calculate. 
The ratio A/gi however is independent of the nuclear species and depends only on the 
position of the nucleus in the molecule. 



B. Degeneracy of arbitrary multiplicity 

The example of the Kramers doublet showed how the paramagnetic shielding can be 
obtained from EPR g- and /L-tensors. In this section we extend this treatment to degenerate 
ground states of any multiplicity and derive a formula for cr p which is the generalization of 
Eq. (ED. 10 

The starting point is again Eq. ffl2l . We have seen that the matrices of m 8 and Ti in a 
Kramers doublet can be reproduced by operator equivalents, linear in the spin components, 
working in a fictitious spin doublet. An expansion of this kind is an application of a general 
theorem to the effect that any n x n matrix can be reproduced by a unique spin operator, 
polynomial in the spin components S^.— Mathematical details about the construction and 
properties of this spin operator are given in Appendix D The main result is that the matrix 
of an operator X in a basis {ipi}i=i...n can be reproduced by a spin operator X$ in a fictitious 
spin multiplet {\S M)} whose dimension equals n, i.e., 2S + 1 = n. Specifically, 16 

25 k 
k=0 q=-k 

where are complex coefficients, unique for given operator X and basis functions ipi. 
The £f ) are irreducible tensor operators, which are basically products of k spin components 
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Si adapted to spherical symmetry.— The formula for the paramagnetic shielding, Eq. (|T2|) . 
depends on six matrices, one for each cartesian component of the magnetic moment m and 
the hyperfme field J- . We now apply the spin-decomposition to each of these matrices and 
write the corresponding spin operators 

2S k 
k=0 q=-k 



,s , ( 15 ) 



k=0 q=—k 

where the Q-numbers have been renamed g^ 1 for the magnetic moment and for the 
hyperfine field, in accordance with the usual notation in EPR spectroscopy. Note that in 
the present formalism the "g-factors" and "A-factors" are complex numbers whose complex 
conjugates are given by 

These relations make the operators in Eq. ( jT5l) Hermitian (See Appendix, Eq. ( 1A.8I) ). The 
spin Hamiltonian 

i i 

= ^ £ £(-1)^4?^ + £ £(-i) ff ^M?4 

i kq i kq 

gives a complete description of an EPR experiment in the degenerate manifold. It is the 
generalization of Eq. (|T3|) to a manifold of arbitrary degeneracy. 

We can now proceed to express the paramagnetic shielding tensor in terms of the EPR 
parameters. In Eq. (TT2]) . we replace the matrix elements in real space by the corresponding, 
and identical, matrix elements in fictitious spin space: 

< = - £ (SM\mi,s\SM')(SM'\^ s \SM) 

M,M' 

= ^TrK,^, 5 ) 

_ (j.b P Wr n(3+9') n ( fc )4( fc ') Tr /"oWo(fc') 



The trace simplifies using Eqs. (IA.7[) and < \AA$ and we get 



fc=0 q=—k 
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This expression is the generalization of Eq. (JHJ). Note that, whereas the spin Hamiltonian 
parameters are basis dependent (see Appendix,), this is not so for the contraction of the 
parameters appearing in Eq. (1161) . The tensor cr p is thus an invariant property of the 
degenerate manifold, as expected. Another such invariant is the magnetic susceptibility 
tensor, which is proportional to the same contraction but with g^- substituted for . For 
the twofold degenerate Kramers doublet this particular contraction is in fact well known as 
the G tensor: G = g T g, with g as in Eq. (TT3T) . 15 



An example of the application of Eq. (fTEj) to the cubic quartet electronic ground state of 



10. In 



certain paramagnetic lanthanide ion impurities in CaF 2 crystals can be found in Ref. 
the next section we apply the formula to the icosahedral sextet that can arise in endohedral 
metallofullerene Ln@Ceo- 



1. Example: paramagnetic shielding in the icosahedral sextet 

The highest degeneracy that can be generated by point group symmetry is sixfold and 
corresponds to the W irreducible spin representation of the icosahedral double group /*.— 
A ground state of this symmetry type may arise from the crystal field splitting of the J 
ground multiplet of the rare earths Ce 3+ , Nd 3+ , Sm 3+ , Dy 3+ , Er 3+ , and Yb 3+ , when placed 
in an environment of icosahedral symmetry— We consider the paramagnetic shielding at 
the central nucleus. 

The high symmetry reduces the number of free parameters in the spin Hamiltonian dras- 
tically. To make maximum use of the symmetry it is best to have the cartesian axes coincide 
with the twofold axes of a Di subgroup of /. This way the x, y, and z directions are 
equivalent,— and the shielding tensor is isotropic and diagonal: <x p = cr p l. The fictitious 
spin representation of W has 5 = 5/2. A vector operator equivalent working in this space 
has to transform as the Tj irreducible representation of the icosahedral group.— This means 
that the linear combinations of tensor operators in Eq. f[T5"j) must transform as the cartesian 
component of T^ Symmetry adaptations of spherical harmonics were published by Qiu and 
Ceulemans.— Using their tables one finds that Ti occurs in the subduction of k — 1 and 
k — 5. The spin operators with k — 1 are just the cartesian spin components S x , S y , and S z . 
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The corresponding operators for k = 5 are as follows (Ref.|2l[ Table 5. (3 n)Tn = M±n+my5)) 



rn _ VWi^a ( q(5) q(5)\ Wjp-i^ ( o(5) Q (5)\ yg/^i / q(5) q (5) 
ia:_ lW2"^5 ^-5^1 16^ ^3 ^-3 J —ft- \!>1 3-1 

T - ■• Vfy-i.s f c(5) , o(5)\ .- aVfft,! / Q (5) Q (5)\ , .. -^g-T.i / q(5) , Q (5) 
^ ~~ * 16^ l, D 5 + 3-5^ ? 16V2 l, D 3 + 3_ 3 j+2 16 + <=>-l 

rn \/35 ( o(5) , «(5)\ , V21 ( c(5) , q(5)\ 3 q(S) 

These expressions are normalized in the sense that the squares of the coefficients sum up 

to one. It remains now to choose a reduced matrix element for S^ 5 \ Eq. ffT6]) adopts a 

simplified form if we choose 

|(g||g( 5 )||g)| 2 _ Kg]gW|]j)P _ (2S+l)S(S + l) _ 35 
11 3 3 2 ' 

for S = 5/2. The second equality follows from our choice of as in Eq. (IA.2j) . 18 

A completely general spin Hamiltonian for Zeeman and hyperfine coupling (with the 

central nucleus) in the icosahedral sextet thus contains just four parameters and is given by 

H s = /i B (<?S + g'T) ■ B + (AS + A'T) ■ I. 

The expression for the paramagnetic shielding follows now readily from Eq. ([16]) : 
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Appendix: Spin operator equivalents 

We consider the following problem. Given an arbitrary, complex square matrix of di- 
mension 2S + 1, to construct a spin operator, polynomial in the components of S, whose 
matrix representation in the fictitious spin manifold \S M) is the same as the given matrix. 
The main results in this appendix are not new and were published by Griffith.— iiZ.>22 For 
convenience, we collect them here in a coherent presentation, and expound in more detail 
on the derivations. 
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1. Orthogonal matrix expansion 



The first non-trivial case is that of a complex 2x2 matrix. This corresponds to fictitious 
spin S = 1/2. It is not hard to see that any such matrix can be expanded as 

a\ + bcr x + ccr y + dcr z , 

where I is the 2x2 unit matrix and the er^ are the Pauli matrices: 

'■'\ioj- CT " = (°~oj' CT "(o-°i)' 

The four complex quantities a, b, c, d are uniquely determined by a set of four linear equa- 
tions, one for each element of the given matrix. The corresponding spin operator is given 
by 

a + 2bS x + 2cS y + 2dS z . (A.l) 

We have seen an example in section II V Al where matrices of Zeeman and hyperfine cou- 
pling in a Kramers doublet were represented by a spin Hamiltonian (Eq. (TL3]) ). which is an 
application of Eq. flA.ll) . (with a = 0, as the matrices considered there are traceless). 

When we move to matrices of dimension n x n (n > 3), the operator equivalent in Eq. 
(lA.ip is no longer adequate to reproduce an arbitrary matrix. This is most easily seen from 
the requirement that the number of parameters in the spin operator must be at least n 2 (i.e., 
the number of elements in the arbitrary matrix), whereas there are only four of them in Eq. 
(lA.lj) . The most obvious way to increase the number of parameters is to include terms of 
quadratic, cubic, and higher degrees in the spin components Si, until the required number of 
parameters is reached. This procedure suffers from a drawback though, which is that linear 
dependencies arise among the terms. For example, there are nine quadratic terms: S x , S x S y , 
S y S x , . . . , but four linear combinations of them are already contained among the terms of 
lower degree (i.e., in Eq. (1A.1I) ). These are S x + Sy + S 2 = S 2 = + 1), S x S y — S y S x = iS z , 
SyS z S z Sy = %S X , and S Z S X S X S Z = iSy. The five remaining linear combinations are 
linearly independent from terms of lower degree and can be added to Eq. (lA.ip to obtain a 
spin operator with nine parameters in total, which allows to expand any 3x3 matrix. 

This procedure may be repeated for higher degrees and one will always find that, among 
the 3 k spin products of degree k, only 2k + 1 are independent from the products of lower 
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degree. These 2k + 1 operators form an irreducible tensorial set of rank k. To proof this we 
make use of the general theory of irreducible tensor operators (ITO) and angular momentum, 
as can be found for example in Zare.— 

Let S q (q = —k, . . . , k) denote the irreducible tensor operators of rank k constructed 
from the products of k spin components S{. is simply a constant, which may be taken 
to be 1. The three components of are (up to a common, real factor of choice, see below) 

So* 1 = S z , S±l = -j={=fS x — iS y ). (A. 2) 

The tensors of higher rank are obtained as follows. is found by coupling with itself 
into a rank-2 ITO. Next, is found by coupling with into a rank-3 ITO, and so 
on. This process may be repeated several times to obtain the ITO of desired rank. We thus 
define our Si k) by Eq. (EO]l and the following recursion relation. 

Sf = N k (-iyJ2 ( k " 1 1 ] Sg-VsV (A.3) 



Qi <?2 -q, 



The multiplicative constant Nk may be chosen freely for each k but we require it to be real, 
in order for the tensor operator to obey certain simple relations under Hermitian and time 
conjugation, to be explained in Section l"31 below. The large bracket is a 3j symbol.— 

We write for the matrix representation of Sq in the basis \SM). The Wigner- 
Eckart theorem tells us that the matrix elements factor into a 3j symbol and a reduced 
matrix element as follows: 



{SM\S^\SM') = (-l) s - M (S\\S w \\S) 



S k S 
-M q M' 



Note that the triangle condition implies that the matrices are zero for all k > 2S. We can now 
show that the (k = 0, 1, ... , 2S) form a basis for the vector space C nxn [n = 2S+1) of all 
complex nxn matrices. First, observe that the set consists of ^? =0 (2A;+1) = (25* + 1) 2 = n 2 
matrices, in agreement with the dimension of C nxn . Second, the matrices constitute a 
linearly independent set, which can be shown as follows. Define an inner product between 
two matrices by imagining each matrix to be reorganized as a vector by sticking the rows 

2 

together and then taking the usual scalar product in C n : 



:A,B) = ^A^ = Tr(AtB). 
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We have now 

Sf\sJ°) = Tr (sf^Sj ) = Y,( S M\S^\S M')*(S M^'^S M') 



MM' 



IS S k\ I S S k'\ 
mm' \M' -M q) \M> -M q'J 

= 2fe + 1 <W<V> (A.4) 

which shows that the Sj^ are orthogonal and therefore linearly independent. Note that the 
derivation leading to Eq. (1A.4|) makes nowhere use of the specific form one may choose for 
the S g k \ such as the one in Eq. (1A.3I) for example. Rather, the result follows straight from 
the Wigner-Eckart theorem, which is by definition true for all ITO's. 

It follows that any given n x n matrix X has a unique expansion in spin matrices of 
the ITOs S { q k) in the basis \S M) (2S + 1 = n). The expansion coefficients are obtained by 
orthogonal projection, using Eq. (1A.4|) . 

2S k 

X = EE (-^Q { q k) S {k) q (A.5) 

fc=0 q=-k 

Q?) Tr(sLfx), a(k) 2k + 1 



\(s\\sw\m 



2. The choice of basis 



We have established that to every square matrix X there corresponds a unique™ set of 



coefficients such that Eq. (I A . 5 [) is true. These coefficients define the spin operator which 



reproduces X in the fictitious spin manifold: 

25 k 

^^EEHWW. (A.6) 

k=0 q=-k 

However, when X is the matrix representation of a true operator X, the spin equivalent X$ 
representing X is not at all unique. Given an orthonormal basis (^i, . . . ,ip n ) and 

X l3 = (^\X\^), 

any unitary transformation of this basis results in an equally valid representation of X, given 
by a transformed matrix X — > U^XU. The transformed matrix may in turn be expanded 
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according to Eq. flA.5j) . yielding a new set of coefficients Q q and hence a new operator 
equivalent X$- Clearly, an infinite number of different spin operators can thus be generated 
and they all represent the same physical situation. Choosing a basis can be seen as estab- 
lishing an association between the true states and the fictitious spin states: ijj(M) -h- \S M). 
A different association leads to a different spin operator, and a suitable choice often leads 
to a simplified spin operator.— 



3. Implications of Hermiticity and time reversal 

The matrices we wish to represent by spin operators are usually Hermitian and they 
often have a symmetry related to time reversal as well. These two properties will naturally 
impose some restrictions on the spin operator equivalent. Under Hermitian conjugation (f) 
and Kramers's time reversal operation the Sq transform as follows: 



5 (fc)t = 

S (k)* = (-ly-kgy 



(A.7) 



which can be verified from Eqs. (1A.2|) and (1A.3I) . using induction on k. 



a. Hermiticity. Suppose now that X is an Hermitian matrix, e.g., the matrix repre- 
sentation of an Hermitian operator in an orthogonal basis of wave functions. There exists a 
unique operator X$, of the form flA.6j) . whose matrix in \S M) equals X. The Hermiticity 



of X then implies 

(SM'\X S \SM) = (SM\X S \SM')* = (S M'\X* S \S M), 

for all M, M', which can only be true if X' s = Xg- Thus the Hermiticity of an operator 
in true space translates into the Hermiticity of the equivalent spin operator, and this is 
independent of the choice of basis in true space. The fact that an Hermitian true operator 
must be represented by an Hermitian spin operator may seem obvious, but we mention 
it explicitly because the behavior of the true operator under time conjugation does not 
automatically lead to the same behavior of the spin operator, as we shall see below. From 
Eqs. (1A.6j) and ( 1A.7j) follow the conditions on the coefficients Q q k ^ that make Xg Hermitian: 



Xl = X s # Q q k >* = (-l)V { -l (A.8) 
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b. Time reversal. The operation of time reversal works on states and on operators. 
We consider the implications hereof on the spin operator equivalent. Let X be the operator 
working in a set of n states. We assume that this set is closed under time reversal. This 
means that there exists a unitary matrix A such that 

n 

(a.9) 

If ijji describes a system of n e electrons then we know that (ip*)* = (—l) ne ipi, and therefore 
AA* = (-l) Ue . From this and the fact that A is unitary follows A = (-l) ne A T , that is, A is 
symmetric for even and antisymmetric for odd number of electrons. In spin space, the time 
reversal operator is given by K$ = Y$Kq, where Y$ is the operator for a rotation through ir 
about the y axis and K takes the complex conjugate of numbers, but leaves the basis kets 
\S M) unchanged.— 

\SM)* = K\SM) = Y S \SM) = {-1) S - M \S -M). (A.10) 

We shall consider the case where the operator X has a definite parity under time reversal: 
X* = eX, with e = ±1. Using Eq. flA.9j) we have for the matrix elements of X 

X*. = Q>i\X\1>jY = mx*\^*) =e^A* ki X kl A lj , (A.ll) 

k,l 

or, in matrix notation, 

X* = eA f XA 

On the other hand, for the matrix elements of Xs we find, using Eq. ( 1A.10I) 

{SM'\X S \SM}* = (S M'\YgXgY s \S M) (A.12) 

We now define As as the spin operator whose matrix representation in \S M) is A. Then we 
can derive from Eqs. flA.111) and flA.12j) and the requirement that Xs reproduces the matrix 
X that 

(S M'\Y^X S Y S \S M) = e(S M'\A S X S A S \S M), 

and therefore 

X* s = eY s A\x s A s Yl (A.13) 

This result shows that, in general, the fact that the true operator X has a definite time 
parity e does not automatically imply that the same is true for the fictitious spin operator 
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Xs- In particular, the behavior of Xs under time conjugation is seen to depend, through 
As, on the choice of basis in the manifold of true states. 

Is it nevertheless possible, by a suitable choice of basis, to obtain that X s = eX^? This 
would require, from Eq. (IA. 13[) . that As = Y$. Then we would have 

(— l)" e = AA* = YsY*s = Y s = (-1) 25 , (A.14) 

where the third equality follows from the reality of the matrix elements of Ys, see Eq. 
(1A.10I) . and the last equality follows from the fact that Y$ is a rotation through an angle 
of 2n, which multiplies a spin ket by +1 if S is integer and by —1 if S is half-integer. We 
conclude from Eq. I I A. 141) that in cases where n e + 2S is odd, we cannot in general obtain 
that X* s = eX s . 

Let us first then consider the other case, n e + 2S even. Recall that the fictitious spin 
quantum number S is determined by the degeneracy of the manifold n = 2S + 1. Thus, 
n e + 25 is even for odd degeneracy in an even electron system and for even degeneracy in 
an odd electron system. For these cases it is always possible to find a basis ip(M) in true 
space satisfying ip(M)* = (—l) s ~ M ip(—M). With this choice we have A = Ys and thus a 
spin operator which has the same time parity as the true operator: X s = eXs- If Xs is also 
Hermitian we find, using Eqs. (1A.6j) . (1A.7I) and (1A.8j) . 



, , N e = +1 and k odd 

X\ = X s and X* s = eX s Q { q k) = for <^ (A.15) 

I e = — 1 and k even 

This result is equivalent to the well-known rule that a spin Hamiltonian can contain only 
even powers of the spin if it represents a time-even operator and only odd powers of the spin 
if it represents a time-odd operator. However, this rule is only valid for a particular choice 
of basis, namely one which transforms under time reversal in exactly the same way as the 
associated spin basis. 

It remains now to discuss the case n e + 2S odd. Although this case is far less common 
than the previous one, it is worth including because the rule mentioned in the previous 
paragraph does not apply here. First note that we may limit the discussion to even electron 
systems having even degeneracy because odd electron systems always belong to the previous 
case, by Kramers's theorem. We have already proven that it will never be possible to have 
the spin operator satisfy X s = eXs, no matter which basis is chosen. Nevertheless, as shown 
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by Griffith,— a simplification of the spin operator similar to Eq. ( 1A.15j) obtains if the basis 
consists of time conjugate pairs. Given a real basis ipi (ip* = ipi) (and this is always possible 
for even-electron systems), construct a new basis ip(M) as follows: 

V(±S) = -^i±# 2 ), 

^(±(5-l)) = i=(^3±#4), 

and so on. Then evidently ip(M)* = ip(—M), and it is easily verified that this corresponds 
to As = (— l)~ Sz ~ s Ys- From Eq. flA. 13[) follows the time-reversal behavior of X$: 

X* = e (-l) s *X s (-l)- s * (A.16) 

Formally, X$ undergoes a rotation through tt about the z axis followed by multiplication with 
e. To translate Eq. (IA.16I) into a rule on the coefficients we use that (— l) Sz Sq k \— l)~ Sz = 
(-l) q S ( q k) and find 22 



Xl = X s and X* = e(-l) s 'X s (-l)- s ' & 



Qj fc) = for 



e = +1 and k — q odd 
e = — 1 and k — q even 
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